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In our previous publication [Kogan et al, Phys. Rev. 48, 
9404 (1993)] we considered the issue of statistics of radiation 
diffusively propagating in a disordered medium. The consid- 
eration was in the framework of diagrammatic techniques and 
a new representation for the intensity distribution function in 
terms of connected diagrams only was proposed. Here we use 
similar approach to treat the issue of statistics in the regime 
of the crossover between ballistic and diffusive transport. We 
find that even small contribution from coherent component 
decreases by one half the intensity distribution function for 
small values of intensity and also produces oscillations of the 
distribution function. We also apply this method to study 
statistics of fluctuations of wave functions of chaotic electrons 
in a quantum dot in an arbitrary magnetic field, by calculat- 
ing the single state local density in the regime of the crossover 
between the orthogonal and unitary ensemble. 

PACS numbers: 42.25-p, 72.10-d, 73.20.Dx, 78.20.Dj 



I. INTRODUCTION 

When electromagnetic wave propagates in a random 
medium and undergoes multiple scattering the density 
of energy in a given point (intensity) / =< E'^ > is 
a strongly fluctuating quantity. In a diffusive regime 
these fluctuations are traditionally described by Rayleigh 
statisticad. The distribution function Pr{I) is a simple 
negative exponential: 
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which corresponds to the equation for the moments: 

< /" >^nl<I >" . (2) 

In Eqs. (0) and (||) the averaging is with respect to 
all macroscopically equivalent but microscopically differ- 
ent scatterers configurations. By diffusive we mean the 
regime when it is possible to discard the " coherent" com- 
ponent of the field, that is 



< 0. 
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This assumption is true for the case r/£ ^ 1, where r is 
a distance from the source to the observation point, and 
^ is a mean free path. The aim of the present publication 
is to get the distribution function for arbitrary relation 



between r and i, when generally speaking the "coherent" 
component of the field should also be taken into account. 

Let us remind how the issue of intensity statistics is 
treated in|-t|he framework of the traditional diagrammatic 
techniqueoB. For the sake of definitness consider the case 
when the point source is situated at the origin; the in- 
tensity is measured in the point r. In the diagrammatic 
representation < / > is given by the diagrams with a pair 
of wave propagators Gq^ and G^^, summed with respect 
to all possible interactions with the scatterers. The n-th 
moment < /" > is given by the set of diagrams with 
n propagators Gq:^ and n propagators Gq^. For us is 
important the following property of diagrams: if the di- 
agram consists from several disconnected parts then the 
contribution of that diagram is equal to the product of 
contributions of disconnected parts, so all the moments 
can be expressed at least in principle through the con- 
tributions of connected diagrams only. To get Rayleigh 
statistics one should ignore all connected diagrams save 
those consistiHg from a pair of propagators (advanced 
and retarded|)3l3. 

Previousljcl we have formulated a perturbation the- 
ory which systematically takes all connected diagrams 
into account. It is connected diagrams with more than 
two propagators which give the deviation from the sim- 
ple Rayleigh statistics for large values of intensity (for the 
tail of the distribution function). In this paper we would 
be interested not in the tail but in the " body" of the dis- 
tribution function. (Because the deviation appears for 
the values of intensity which are inversely proportional 
to r this limitation looks quite natural in the regime con- 
sidered). Hence we would discard connected diagrams 
with more than two propagators. On the other hand pre- 
viously we ignored the diagrams which contain isolated 
(dressed) propagators. The ground for this is that such 
diagrams give contribution of the order of exp{—r/£), and 
in the fully developed diffusive regime one can ignore 
them. Here we should generally speaking take them into 
account. 



II. CROSSOVER BETWEEN DIFFUSIVE AND 
BALLISTIC REGIME 

For < > we consider only diagrams which consists 
of isolated (dressed) propagators and pairs of propagators 
(advanced and retarded one). The sum of all diagrams 
can be written down in the following way: 
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where A = \< E >\\ B =< I > -\ < E > , and 
P{n,'m) is the coefficient of purely combinatorial origin: 

(n!)2 



P{n, m) = 
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The distribution function P{I) is connected with its 
moments by usual way: 
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Substituting the Eq. (Q) into the Eq. (g) and changing 
the order of summation,which is possible due to absolute 
convergence of the series, we get: 
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Taking into account that for n > m 



(n-m)! 

we finally obtain the distribution function in the form 



(8) 
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(the intensity is measured in the units of < / > and 
parameter a = A/{A + B) ^ \ < E > / < I >, is 
introduced) . It is obvious that for a = (purely diffusive 
transport), the Eq. (^) coincides with Eq. (|l|), so this 
parameter gives the degree of deviation from Rayleigh 
statistics. In the opposite limiting case a = I Eq. (g) 
gives P {!) — S{I—1). This limiting case means that 
there are no fluctuations, which can occur only when 
there is no scattering at all, in other words in the ballistic 
regime. 

For intermediate value of a the distribution function 
can easily be calculated numerically. This behavior is 
represented on Fig. 1. We see that even for a <C 1 
(almost diffusive transport) taking into account of the 
" coherent" component of the field drastically changes the 
distribution function for small /; in particular for / = 
we get: P = 1/2. Also interesting are the predicted 
oscillations of the distribution function. 



III. CROSSOVER BETWEEN ORTHOGONAL 
AND UNITARY ENSEMBLE 

Now we want to consider a second problem, which 
though physically different from the considered above is 



as we would see very similar mathematically. It's the 
problem of statistics of fluctuations of wave functions of 
chaotic electrons in a quantum dot in an arbitrary mag- 
netic fi^ld- This problem has been a subject of several 
studiesaiO. The distribution function P^ (I) of a single 
state local electron density / = (?") P was shown to 
be for unitary ensemble (strong magnetic field) a simple 
negative exponent (like in Eq. (|^)), and for orthogonal 
ensemble (no magnetic field) 
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In the papei0 the distribution function was calculated 
for the case of arbitrary magnetic field using the super- 
symmetry technique. Let us see what the simple method 
gives in application to this problem. Both limiting cases 
can be easily obtained using the method of moments, if 
we postulate for the density moment 



< /" >=< ^ . . . ^ tp* ...tp* > 



the existence of Wick theorem. Then we obtain 
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where the multiplier nl is simply the number of possible 
pairings between amplitudes, and hence Rayleigh statis- 
tics For orthogonal ensemble the wavefunction is real, 
that is there is no difference between tp and V'* a ny two 
propagators can be coupled and instead of Eq. ( p^ we 
get : 



<V' •■•'(/'>= (2n-l)!! <i:^P 
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which gives us Eq. (|T^). (This statistics, by the way, 
is also known in the theory of propagation of classical 
waves in a disordered medium; Eq. (|l^) gives distribu- 
tion function for the case of acoustic waves). This gives 
us the idea how the general issue of electron density dis- 
tribution function should be addressed. We should take 
into account both the averages < tpip* > and < i/jip > 
(not supposing of course them to be equal). Then tak- 
ing into account all possible ways of coupling and using 
simple combinatorics we get: 

[n/2] 
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where P(n, m) is the coefficient of purely combinatorial 
origin: 
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Substituting Eq. (Q into Eq. (g) and changing the 
order of summation we get: 
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Using Eq. (^) we obtain: 

where the density is measured in the units of < / > and 
Y =< tpip* > / \< tpi/j >| — 1 is the parameter which 
shows the character of the ensemble (for orthogonal en- 
semble Y = and for unitary Y = oo). It is obvious that 
for y = oo we get Eq. (^ and for F = we get Eq. (10) 



For intermediate value of Y the distribution function 
can easily be calculated numerically. To compare our 
result with the results of Ref. on Fig. 2 we plotted the 
distribution ip (t) = 2rPe (t^ as function of X = ^/Y 
(looking at / ^ F ^ asimptotics of Eq. (^7|) we 
see that apart from numerical coefficient very close to 
1 our parameter X is equal to that introduced Ref. |^). 
Comparing our Fig. 2 with Fig. 1 of Ref. we see that 
they look very much alike, though analytical expression 
for distribution function we have got differs from that of 
Ref. 0. 



ACKNOWLEDGEMENTS 



The authors acknowledge the financial support of the 
Israeli Academy of Sciences. 



^ J. W. Goodman, Statistical Optics (J. Wiley, New York, NY, 
1985). 

^ B.Shapiro, Phys. Rev. Lett, 57, 2168 (1986). 
^ E. Kogan, M. Kaveh, R. Baumgartner, and R. Berkovits, 
Phys. Rev. B 48, 9404 (1993); Physica A 200, 469 (1993). 
■^F. Wegner, Z. Phys. 36, 207 (1980). 

^ B. L. Altshuler and V. N. Prigodin, Zh. Eksp. Teor. Pliys. 

95, 348 (1989) [Sov. Phys. JETP 68, 198 (1989)]. 
^ K. B. Efetov and V. N. Prigodin, Phys. Rev. Lett, 70, 1315 

(1993) . 

^ V. I. Fal'ko and K. B. Efetov, Phys. Rev. B 50, 11 267 

(1994) . 



FIG. 1. Distribution function for: (1) a = .1 (solid line), 
(2) a — .5 (dashed line), (3) a — .9 (dot-dashed line). 



FIG. 2. The distribution function <p{t) of the local ampli- 
tude T shown in the crossover regime for different values of 
the parameter X. 
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